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Abstract 

We discuss the non-Boltzmannian nature of quasi-stationary states in the Hamil- 
tonian Mean Field (HMF) model, a paradigmatic model for long-range interacting 
classical many-body systems. We present a theorem excluding the Boltzmann-Gibbs 
exponential weight in Gibbs T-space of microscopic configurations, and comment a 
paper recently published by Baldovin and Orlandini (2006). On the basis of the 
points here discussed, the ongoing debate on the possible application, within ap- 
propriate limits, of the generalized ^-statistics to long-range Hamiltonian systems 
remains open. 
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1 Introduction 



The Hamiltonian Mean Field (HMF) model is a paradigmatic toy model for 
long-range interacting Hamiltonian systems [TJ. In a recent paper [2], Baldovin 
and Orlandini presented new molecular dynamics numerical results for the 
HMF model in contact with a thermal bath. They conclude that the energy 
distribution in its quasi- stationary state (QSS) is of the Boltzmann-Gibbs 
(BG) form. In view of their own numerical results, we cannot agree with their 
conclusions. In this paper we present three points that are relevant for the 
interesting questions raised in [2]. 



2 Discussion 

In what follows we present three important points regarding the non-Boltz- 
mannian nature of the QSS in the HMF model. 

(i) The nature of the QSS in the canonical setup presents differences with that 
in the microcanonical setup. 

It has already been demonstrated that, in the microcanonical ensemble, the 
anomalous dynamics (concerning the presence of a glassy phase, hierarchical 
structures, velocity correlations, aging, vanishing Lyapunov exponents, etc), 
and therefore the statistical description of the QSS, is strongly dependent on 
the initial conditions PPf5] . A fortiori, it is thus not at all justified to extrapo- 
late the conclusions for the canonical QSS studied in [2116] to the microcanonical 
ones. Although the authors themselves state in [6] that the canonical QSS's 
are only reminiscent of the microcanonical ones (in fact they differ in the life- 
times, the relaxation times, and, most probably, also in the correlations), they 
neglect it in [2]. In order to stress this point, in Fig. 1 (a) we plot three numer- 
ical simulations in the microcanonical ensemble (averaged over 10 events) for 
U = 0.69, iV = 5000 and for different initial magnetizations Mq. Water-bag 
initial conditions were used |3j . Following the Central Limit Theorem aspects 
recently revealed in [7], we illustrate in Fig. 1 (b) the non-Boltzmannian na- 
ture of this QSS regime by studying the probability density functions (pdfs) 
of the stochastic values y defined as follows: yj = Yh=i f° r 3 = 1' 2, L, 
where lj (i) are the momenta of the ji/i-rotator taken at fixed intervals of time 
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Fig. 1. (a) The temperature as a function of time for the HMF model at 
U = 0.69 and N = 5000 obtained for three initial magnetization values, respec- 
tively Mq = 1, 0.8, 0, and averaged over 10 events. The system remains trapped in 
a QSS at a temperature value lower than the equilibrium one. The latter is drawn 
together with the limit temperature at this energy [3]. (b) The pdf of stochastic 
values y obtained by summing rotators momenta in the QSS regime illustrated in 
panel (a), see text. Different initial magnetizations Mq are reported. A Gaussian 
(dashed curve) and a g-Gaussian with q = 1.75 (full curve) are also reported for 
comparison, (c) For the initial magnetization Mq = 0.8, we plot the pdfs obtained 
changing the time interval for recording the momenta of the rotators. No significant 
change is observed. See text for further details. 
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i inside the QSS regime. We considered L = 5000 and n=500 intervals, one 
every 10 units of time. The length of the time interval was varied as shown in 
Fig. 1 (c) for the initial magnetization M = 0.8 without any relevant change 
in the pdfs. An average over 10 events was also considered, in order to have 
50000 values for the y variables for each pdf. In case of no global correlations, 
as one would expect at equilibrium, the resulting pdf should be a Gaussian. 
As Fig. 1(b) shows, in our case, the Gaussian distribution does not occur. In 
fact we have pdfs whose shape depends on the initial magnetization M and 
the central region (up to a few mean standard deviations) of most of them 
is well fitted by g-Gaussians. In particular for M = 0.8 we get a g-Gaussian 
with q = 1.75 as in [7j. We do not think that this is a universal value in this 
case, since we do not have a strictly zero largest Lyapunov exponent as in 
[7]. For the HMF model, the latter goes to zero as 1/iV 1 / 9 in the QSS regime 
[3], while it remains finite at equilibrium at this energy value. We normalized 
the data to the corresponding standard deviations in order to compare them. 
This comparison clearly shows that strong correlations are present in our case 
indicating a deviation from the standard Central Limit Theorem, as in [7]. A 
more complete study on this point will be reported elsewhere. In any case the 
arguments discussed in [2], and in particular in Fig. 1(a) of [2], do not apply 
to the microcanonical QSS. The effects become increasingly stronger when the 
initial magnetization gets closer and closer to 1. 

(ii) The energy distribution of the QSS cannot be of the Boltzmann-Gibbs ex- 
ponential form. 

As clearly seen in Fig. 1(a) of [2J, the one-body momentum distribution is non- 
Gaussian. This excludes the BG exponential form as the energy distribution 
in full phase space T, as proved (by reduction to the absurd) in what follows. 
This analytical fact can easily be enlarged to n-dimensional rotators, n > 1: 
n = 1, 2, d, and oo correspond respectively to the Ising, XY, Heisenberg, and 
spherical models. In what follows we illustrate it, however, for n — 2, the case 
addressed in [2]. 

Theorem 

Let H = K+V be the Hamiltonian of a M-body classical system, where the ki- 
netic energy is K = Y^Lx 21 > ^ > ®> anc ^ ^ ne potential energy V contains arbi- 
trary integrable one-body, two-body, three-body... terms, concerning (isotropic 
or anisotropic) rotators localized at a (irregular or regular <i-dimensional) lat- 
tice (each term might be characterized by a distance-dependent coupling con- 
stant, which can be summable or not in the M — * 00 limit). And let the 
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one-momentum marginal probability distribution 

p(h) = J dl 2 dl 3 ...dl M dQ\ d9 2 ..A6 M p{l 1} l 2 , Im, #2, • @m) 
associated with the joint distribution 

p(h, h, hi, 9i, #2, Om) 

be non-Gaussian. Then 

p(h, h, hi, &i, 02, Om) 
cannot be proportional to 

e ~/3 H(Ii,12,—,Im,8i,02,—,8m) ^ > Q 



Proof: 

Assume that 



p(h, h, Im, @i, @2, •••) 6m) 



Je -,H > 

then 

/ dl 2 dl 3 ...dl M dd\dQ 2 ..A6m e~^ H 
P ^ = j dh dl 2 dl 3 ..Al M dQ x dd 2 ..A9 M e~ pH = 

J dl 2 dl 3 ..Al M e~P K J d6 1 d6 2 d6 3 ..A6 M e~ pv e ~^ /2I ^J di ie ~^^ 



f dh dl 2 dl 3 ..Al M e~P K J d6 l d0 2 d0 3 ..A6 M e-? v " nf£i / dh e-W 21 

e" 



-1311/21 



Jdhe-^i/ 21 ' 

which is a Gaussian, thus contradicting the hypothesis. Q.E.D. 

Hence, the Central Limit Theorem verification exhibited in j2] refers to a nec- 
essary but insufficient property. 



(Hi) The Baldovin-Orlandini results do not exclude the nonextensive statistical 
mechanics q- exponential form. 



Contrary to what is stated in [2], g-statistics [8] cannot be excluded on the 
basis of what they exhibit in their Fig. 3(d). Indeed, in Fig. 2 we plot the 
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Fig. 2. The statistical weights corresponding to BG statistics (full circles), and to a 
typical illustration of nonextensive statistics {q = 1.5; solid line), for the M = 1000 
example exhibited in Fig. 3(d) of [2], see text for details. 

statistical weights corresponding to BG statistics (full circles), and to a typical 
illustration of nonextensive statistics (q = 1.5; solid line), for the M = 1000 
example exhibited in Fig. 3(d) of [2]. We used nearly the same number of 
points of Fig. 3(d) in the interval Ehmf £ [664, 720], which covers 4% around 
its central value 692. The linear correlation is r = —0.99996 (r = —0.99997 in 
[2 J for a similar analysis). The remarkable closeness of the q — 1 and q — 1.5 
examples comes from two facts. The first one is that 



Notice the absence of the term linear in x in the latter expression. While the 
second one is that the spreading around the central value 692 is very small. 
Notice also that we rescaled T = 0.38 by a factor M, which corresponds 
to multiplying the Hamiltonian (1) of [2 J by M and then rescaling time in 
the kinetic energy (see [9]). This makes the microscopic two-body interaction 
of the model independent from the macroscopic quantity M, as desirable. 
Consider also that, in the ^-statistical weight of Fig. 2, we have subtracted 
from Ehmf a characteristic macroscopic energy, namely the central value 692. 
This procedure is long known to be mandatory (see [TO]) for any q ^ 1 in order 
to preserve the invariance of probabilities with regard to an uniform shift of 
the energy scale, i.e., in order to preserve the freedom of the choice of the zero 
of energies. Such subtraction, although allowed for all values of q, is clearly not 
mandatory for the BG exponential weight (q = 1), since the desired invariance 




[1 + (1 - q)x] 



V(l-9) 



e x [l + -(q - l)x 2 + ...], for x -> 0. 
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is anyhow guaranteed by the fact that the exponential of a sum factorizes into 
the product of exponentials (whereas the g-exponential of a sum does not 
factorize into the product of g-exponentials if q ^ 1). 



3 Conclusions 

On the basis of the points here discussed, the results presented in [2] should be 
considered with extreme care, and by no means they close the ongoing debate 
on the possible application, in appropriate limits, of generalized g-statistics 
to long-range Hamiltonian systems. Finally, let us mention for completeness 
that neither the arguments recently advanced in [11] exclude the possible 
application of (/-statistics, especially for initial magnetization close to 1, but 
this point is beyond the scope of the present paper and it will be discussed 
elsewhere. 
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